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Abstract 

We discussed neutrino masses and mixings in SUSY SO(IO) model where quarks 
and leptons have Yukawa couplings to at least two 10 and one 126 Higgs scalars. 
In this model, the Dirac and the right-handed Majorana mass terms are expressed 
by linear combinations of quark and charged lepton mass matrices, which then 
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there are various solutions to reproduce a large mixing angle for — Uj- and a small 
mixing angle for Vf. — as well as the hierarchical mass spectrum of neutrinos. 


* e-mail: oda@phys. wani .osaka-u .ac.jp 
le-mail: takasugi@phys.wani.osaka-u.ac.jp 
le-mail: tanaka@phys.wani.osaka-u.ac.jp 
§e-mail: masaki@phys.wani.osaka-u.ac.jp 



1 Introduction 


Super-Kamiokande group announced the remarkable report [1], the evidence of the neu¬ 
trino oscillation and the neutrino masses based on the atmospheric neutrino observa¬ 
tion. The — Vf. oscillation scenario is excluded by the CHOOZ data[2] and also the 
Super-Kamiokande data[l], and the — Vr oscillation is favored, although the interpre¬ 
tation of the — ^sterile K possible. Within three neutrino scenario, they showed that 
sin^ > 0.8 and Am 23 range 10 ^eV^ ~ 10 ^eV^. On the other hand, from 

the observation of the day and night difference of the solar neutrino flux, it seems that 
the small mixing for z/g — z/^ is favored for the solar neutrino problem[3]. 

In this paper, we consider these facts seriously and seek the scenario to reproduce the 
pattern of the neutrino mixings and the neutrino mass spectrum. In particular, we aim at 
how the large mixing between z/^ and Vr arises, while one keeps the small mixing between 
Vf. and z/^ in the framework of SUSY SO(IO) GUT models. We consider the model where 
fermions have Yukawa couplings to at least two 10 and one 126 Higgs scalars. In this 
scenario, the Dirac and the right-handed Majorana mass terms are expressed by linear 
combinations of quark and charged lepton mass matrices and thus the neutrino mass 
matrix arising from the see-saw mechanism is also determined by quark and charged 
lepton mass matrices. In the basis where the u-type quark mass matrix is diagonal, 
the d-type and also the l-type(charged lepton) mass matrices are expected to be almost 
diagonal so that it is a nontrivial problem how to obtain the non-hierarchical neutrino 
mass matrix for the part related to the second and the third generations by using these 
hierarchical mass matrices, which is needed to obtain the large mixing between z/^ and 

Vr- 

The model which we consider have been discussed intensively to get the unihed de¬ 
scription of quark and lepton masses and mixings. Babu and Mohapatra[4], and Lee and 
Mohapatra[5] considered the minimal (SUSY) SO(IO) GUT model, where quarks and 
leptons have Yukawa couplings to only one 10 and one 126 Higgs scalars, in order to get 
predictions of neutrino masses and mixings. Along this line, the texture zero analysis 
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based on models of one 10 and two or three 126 Higgs scalars was made by Babu and 
Shafi[6], and Achiman and Greiner[7]. All these models have predicted small mixing 
angle for — Vr- The particular interest of these models lies in the fact that neutrino 
mixings and ratios of neutrino masses are predicted. 

Recently, Brahmachari and Mohapatra[8] discussed that minimal SUSY SO(IO) mod¬ 
els with one 10 and one 126 is unable to predict a large mixing angle for — Therefore, 
they considered SUSY SO(IO) models where Higgs multiplets are in two 10 and one 126 
representations. By considering type H see-saw mechanism where neutrino mass ma¬ 
trix consists of the left-handed Majorana mass term and the see-saw term, they found 
a solution which predicts a large mixing angle for and a small mixing angle for 

Motivated by the work in Ref.8, we consider SUSY SO(IO) models where quarks and 
leptons have Yukawa couplings to more than two 10 and one 126 Higgs scalars. The 
model is essentially the same as those by Brahmachari and Mohapatra except that we do 
not consider the left-handed Majorana mass term. We assume that the low energy theory 
of these models is the MSSM with two Higgs doublets which are linear combinations of 
the doublets in the lO’s and the 126. We assume that Hu = Q;ii7„(10i) -|- a2Hu{102) + 
0 : 377 ^,( 126 ) and Hd = /5i77rf(10i)-1-/5277^(102)+/53i7rf(126). Quark and lepton masses come 
from the Yukawa couplings, 

Wy = hi^abi^aAHiilO) + /afe^/’a^feA(l26) , (1) 

where tjja is the 16 dimensional fermions with the family index a. The matrices 7j’s and 
/ are 3 by 3 complex symmetric matrices. Quark and lepton mass matrices are given by 

n 

TTu Vujhj -|- Kuf , Mfi ^ ) ^djhj Kdf , 

j 

n 

MyD = Ej Vujhj - 3Kuf , Mi = Y^ Vdjhj - Sudf , 

j 

Mur = VRf , ( 2 ) 

where Vui is the vacuum expectation value of 77„(10i) multiplied by the ratio of Hu in 
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Hui^Oi). Others are defined similarly. They satisfy 




!](<■ + vjj) + Kl + Kj = 246GeV . 


The parameter is the scale of the right-handed Majorana neutrino masses. 
From Eq.(2), we obtain the relation for n > 2 


= M, - r{Md - Ml) , 
Mur = R{Md — Ml) , 


( 3 ) 


( 4 ) 


where r = Ku/K d and R = VR/And- Since we consider that the neutrino mass matrix 
m,y is derived by the see-saw mechanism as m^, = —M^jjM~^M^r,, is essentially 
determined by quark and charged lepton mass matrices. Only other parameters are the 
ratios of vacuum expectation values, r and R. The parameter R is used to control the 
overall normalization of neutrino masses and therefore r is the only adjusting parameter 
to fix neutrino mixings and the ratios of neutrino masses. The models presented above 
are quite tight so that there needs some mechanism which naturally leads to the large 
mixing between and Ur. 

The topics in this paper are as follows: 

(i) We show that the minimal model consisting of one 10 and one 126 Higgs is excluded. 
This is because of the inability to reproduce experimentally observed pattern of mass 
spectrum of d-type quarks and charged leptons simultaneously. Our reason is severer 
than the reason raised by Brahmachari and Mohapatra[8]. They argued that the model 
is unable to predict a large mixing angle between — Ur and thus the model is rejected. 
The outline of our discussion is given in Appendix A. 

(ii) We consider the (type I) see-saw mechanism to obtain the neutrino mass matrix 
in contrast to Brahmachari and Mohapatra[8] who used the type II see-saw mechanism 
where the left-handed Majorana mass term is added to the see-saw term. We simply 
avoid to introduce an extra freedom due to the left-handed Majorana mass term. 

(hi) We give a qualitative argument how to get the less hierarchical structure of neutrino 
mass matrix which is derived through the see-saw mechanism by using the Dirac and the 
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right-handed Majorana mass terms which are given as linear combinations of hierarchical 
u-type, d-type and the charged lepton mass terms. 

(iv) We found many possible ranges of parameter r to lead a small mixing for z/g — 
and a large mixing for — Vj. as well as the hierarchical neutrino mass spectrum. Each 
region corresponds to slightly different mixing angles and neutrino masses. 

The paper is organazed as follows: In Sec.2, we explain the general structure of the 
Dirac mass matrix and the Majorana mass matrix, and then we explaine the strategy 
how to obtain the large mixing between z/^ and z/^. In Sec.3, we present the numerical 
analysis on whether models can reproduce the atomospheric neutrino as well as the solar 
neutrino experimental data. The summary is given in Sec.4. 

2 The unified description of quark and lepton masses 
and mixings 

In this section, we present the qualitative argument how to derive the large mixing 
between z/^ and Vr- 

2.1 The basis of mass matrices 

There is a freedom to choose the basis of mass matrices. Because neutrino mass matrices 
are expressed as linear combinations of quark and charged lepton mass matrices, we can 
transform all mass matrices simultaneously by a unitary matrix such that = U'^MkU 
where k = u,d, I, uD, uR by a unitary matrix U. Since all mass matrices are symmetric 
ones, we choose U such that it diagonalizes 

Mi = , (5) 

where = diag(|mt,|, Imd, \mt\). Next, we introduce a unitary matrix D which diago¬ 
nalizes Mrf, i.e., D^MdD = Da with = diagdm^l, \ms\, \mh\). By using U and D, the 
general form of the CKM matrix is given by = cpiKcpd where is a special form 
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of the CKM matrix and 0^ (i=u,d) are diagonal phase matrices. By using K and phase 
matrices M'^ is expressed by 

M' = U^M,U = U^D*D,D^U = . (6) 

Finally, we transform all matrices by the phase matrix <pu as This is the 

basis which we use in below: 

Mu = Du, Md = K*bdK^ (7) 

and Ml, and where = (plDuCpi = diag(m„, rric, mj) and = (t)*dDd4>\ = 

diag(mrf, mb). In this basis, all quark masses are complex quantities. Now we have 

M,d = Mu-rM , Mur = RM , (8) 

where 

M = Md-Mi. (9) 

The left-handed neutrino mass matrix rrii, is given by the see-saw mechanism as 

m,y = ■ ( 10 ) 

Thus, neutrino mass matrix is essentially determined once quark and charged lepton mass 
matrices are given. Only other parameters in the model are r and R. The parameter 
R determines the overall scale of neutrino masses. Thus, the parameter r is the only 
adjusting parameter to reproduce the desired neutrino mixing angles and neutrino mass 
ratios. 

2.2 Quark and charged lepton masses and CKM mixings at the 
GUT scale 

We use the following quark and charged lepton masses, and the CKM mixing angles at 
GUT scale (2 x lO^^GeV) which were estimated in the minimal supersymmetric standard 
model with tan/? = 10 following Fusaoka and Koide[9]: 
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\mu\ 

= 0.00104 

\mc\ 

= 0.302 

\mt\ 

= 129 


\md\ 

= 0.00133 

\ms\ 

= 0.0265 

\mb\ 

= 1.00 , 

(11) 

\me\ 

= 0.000325 


= 0.0686 

Inirl ■ 

= 1.171 



sin 012 = -0.2205 , sin ^13 = 0.0026 , sin 623 = 0.0318 . (12) 


where fermion masses are defined in unit of GeV. Although these values of parameters 
should have errors, we neglect errors since our purpose is to answer whether the parameter 
ranges of r exist to reproduced the desired neutrino mixings and neutrino spectrum. The 
above values of CKM mixings at GUT scale are given by taking into account of one loop 
contribution by keeping only nit and nib. 

In the numerical analysis, we use values in Eqs.(ll) and (12). For the qualitative 
analysis, we use the Wolfenstein form for the GKM matrix K at the GUT scale 

A AA^e-*^ \ 


K 


-A + |A 2 A 5 

V AA^ - AA^e*'^ 


1 - - (I + 

-AA2 + lAA^ - AA^e*^ 


AA2 


+ 0 (A®), 


l-iA^AV 


(13) 


where A = 0.2205, A = 0.6540 and A = 1.100 by using the mixing angles in Eq.(12). 


2.3 Explicit form of quark mass matrices at GUT scale 


Due to the hierarchy of magnitudes of quark masses, we parametrize Du and Dd as 
follows. 







\ 

Du = mt 

^ctA" 


,Dd = nib 

UX^ 



1 

ij 


[ 

1/ 


where ^’s are quantities of order unity. From Eq.(ll), we have |Gt| = 0.318, |.^ct| = 0.990, 
lUl = 0.563 and \^sb\ = 0.545. 

By using them, one hnds 


Md 


rUb 


— {^db + 


CsbA^ + {—^sb + A^)A^ 


AA2 - e.feAA^ 


V AA^e-*^ - AA 2 - {.feAA^ 1 - A^A^ J 
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/0(A«) 0(\‘) 0(A')\ 
0(A«) 0(A®) 0(A®) 
1>0(A') 0(A®) 0(\‘)l 


( 15 ) 


It is interesting to observe the difference between the mass hierarchy of u-type quarks 
and that of d-type quarks. While ~ O(A^) and \mclrnt\ ~ 0(A"^), \mdlrnb\ ~ 

0(A"^) and \ms/mb\ ~ O(A^). That is, the mass hierarchy of d-quarks is much less sevierer 
than that of u-quarks. Next, we observe that {Md )22 ~ (dTd )23 ~ O(A^). These are crucial 
in the following discussions. 


2.4 The hierarchy in the neutrino mass matrix 



/ 1 

e e'\ 

0 ~ 

—e 

c s 


1 e' 

—s c / 


Firstly, we discuss what kind of neutrino mass matrix is required from the recent data. 
From the data[1],[2],[3], the neutrino mass mixing matrix O is almost hxed aside from 
CP violation phases as 


(16) 


where s = and c = cos 6*^,- with 6 ^r is the mixing angle between and Uj. 

neutrinos, e ~ A^ and e'/e ~ s/(l + c). From this mixing matrix, the expected neutrino 
mass matrix is given by 

/ ri + e'^ se' ce' \ 

rriy = ODiyO^ ~ niu^ se' + r 2 C^ sc , (17) 

\ ce' sc r 2 s‘^ + c^ / 

where ri = ^2 = ~ ±1/10. If we take sin^26'^T- > 0.7 for the 

experimental allowed region, we see that the submatrix relevant to and the Ur should 
have less hierarchical structure than quark mass matrices as 


0.23 ± 0.77r2 0.42 

0.42 0.23r2 ± 0.77 



0.77±0.23r2 0.42 

0.42 0.77r2 ± 0.23 


(18) 
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The above matrices show variations of their components. They corresponds to the angle 

= 6 * 0 , 7J'/4 and 7 r /2 — 6 *o with sin^ 26*o = 0.7. From the above analysis, we observe the 
followings: (1) The 1-2 and 1-3 elements are of order A^, while the 1-1 element is of order 
A"^. (2) In general, the non-hierarchical structure for the part related to and z/,-. Even 
in the extreme case, the hierarchy is at most of order A. 

2.5 The mechanism which leads the large mixing between 
and Ur 

In the basis where is diagonal, takes a hierarchical form. Then, it is natural to 
suppose that Mi also take a hierarchical form. On the other hand, in order to get the 
non-hierarchical neutrino mass matrix for the part related to z/^ and Ur, at least one of 
M^d and M^r should take the non-hierarchical form for the relevant part. Since 
and M^r are linear combinations of hierarchical mass matrices. Mu, Md and Mi, there 
needs some mechanism to get the non-hierarchical structure for M^d and/or M^r for 
the relevant part. This is a necessary condition and does not imply the desired form of 
neutrio mass is obtained. However, we seek this possibility. 

The hint lies in the fact that the 2-2 element of Md is the same size as the 2-3 element. 
We consider how to obtain the non-hierarchical form of M^d- By adjusting r, we make 
the 3-3 element of M^d = Mu — rM as small as of order X^rrit. Then, rric and rUu do 
not contribute to because of the large hierarchy of u-type quark masses. Thus, only 
rrit in contributes to in M,^£,. Thus the non-hierarchical suturucture arises with the 
above condition for the 3-3 element. 

We consider the above condition in detail, which is treated as two separate cases. 

(a) mbnir < 0 case 

In this case, M^r has a hierarchical form as we see from Eqs.( 8 ) and (9). We require 
that My£, has a non-hierarchical form for the relevant part. This is achieved by requiring 
that the 3-3 element of M^q to be of order of the 2-3 element: 


rrit — r{mb — rrir) ~ . 


( 19 ) 



That is, the value of r is hxed to be of order 

r -- 50 . (20) 

mb — mr 

In particular, the range of r which gives a non-hierarchical neutrino mass structure will 
consist of two parts. One is the region r < 50 and the other is the one r > 50, because 
the exact equality mt = r{mb — mr) gives the vanishing 3-3 element of neutrino mass 
matrix so that we can not reproduce the desired mass matrix given in Eq.(18). 

(b) mbmr > 0 case 

In this case, MuR has a non-hierarchical form due to the cancellation between mb and 
mr- There are two cases. 

(b-1) The case where Myo has a non-hierarchical form: 

r - ^ - 750 . (21) 

mb — mr 

The range of r which gives a non-hierarchical neutrino mass structure will consist of two 
parts, the one r < 750 and the other r > 750. 

(b-2) The case where has a hierarchical form: 

mt 

r « - . 

mb — mr 

The above conditions are of course only the necessary ones to achieve a non-hierarchical 
form of the neutrino mass matrix. The problem is whether there are regions of r which 
reproduce the desired neutrino mixing angles and more importantly the hierarchical spec¬ 
trum of neutrino masses. There is no guarantee of the existence of such parameter region 
of r. We have to calculate the neutrino matrix for a given r and examine whether the 
desired mixings and mass ratios are realized, especially by paying a special attention on 
the above mentioned parameter regions of r. 

3 The analysis 

Before the analysis, the following comments are in order: 
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(i) The parameters 

In addition to r and i?, we have sign freedoms of quark and lepton masses. Since one 
of the phase can be hxed, we choose nib > 0. If we scale M^d by mt, the parameter r 
enters as r/rrit. Thus, we can £x r > 0, while we allow rrit to take both positive and 
negative signs. Thus, we fix nib > 0 and r > 0 and take all combinations of signs of other 
fermion masses. 

(ii) The desired neutrino mixings and masses 

We consider the following constraints on ranges of neutrino mixings and neutrino 
masses[l],[10]: 


3 X 10“^ < sin^ < 2.0 x 10"^ , 0.7 < sin^ 20^^ , 
3 X 10-* < ^ < 5 X 10-2 , 


(23) 


Amis 

where we used 3xlO“®eV^ < Anil 2 < lxl0“^eV^ and 2xlO“"^eV^ < Am^g < lxl0“^eV^. 
Since we are dealing with the hierarchical mass spectrum of neutrino case, the overall 
normalization is fixed by the mass ni,yr- The parameter R is determined to £x ni,yr in 
the range of 


3 X lO-^eV^ < Amj^ < 1 x lO'^eV^ , 

2 X lO-^eV^ < Amig < 5 x lO'^^eV^ . (24) 

In the above, we defined Am|^ = ml — and ml « m^ « m\. 

(iii) The mixing angles 

Since we are looking for the solutions which reproduce the hierarchical neutrino mass 
spectrum and the mixings in Eqs.(23) and (24), we may treat the three neutrino mixing 

as if it is due to the two neutrino mixing. In other words, we may define angles by 

sin6*e;^ ~ {Uy)i 2 and sin^^,- ~ {U^) 23 , where Ui, is the neutrino mixing matrix. With this 
approximation, we seek the range of parameter r which reproduce the desired mixings 
and ratios of masses to solve the atmospheric neutrino problem and the solar neutrino 
problem. Once we find the solution, then we examine whether the mixings and masses 
that we obtained do not violate the CHOOZ bound for the Ug — 9x oscillation[2]. Since 
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Am ^2 << ^^23 — ^^13 — ■'^ 3 ) can treat the mixing angle as sin6*ex = (h^iy)i 3 and 
the mass squared difference as Am^ ~ in our three neutrino mixing scenario. 

3.1 The simplification of the problem 

In this paper, we use the following simplihcation; 

(i) The CP violation is neglected. Explicitly, we perform the numerical analysis by setting 
the Kobayashi-Maskawa CP violation angle 6 zero and taking quark and charged lepton 
masses are real. However, we set signs of fermion masses are free so that we have to 
consider all combinations of signs of quark and charged lepton masses. 

(ii) Ml is assumed to be diagonal in the basis we adopted in this paper. 

With these simplihcations, neutrino mass matrix rrii, is determined by quark and 
charged lepton masses, CKM mixing angles, r and R. The parameter R determines the 
overall scale of so that it plays a role of adjusting the overall scale of neutrino masses 
in the see-saw mechanism. The neutrino mixing angles and the ratios of neutrino masses 
are solely determined by only one parameter r. Thus, the present model is quite tightly 
constrained one. 

Once the ranges of parameter r which reproduce the desired neutrino mixing angles 
and ratios of masses are found, the introduction of CP violation will relax the region of 
r and R. The same is expected if we relax (ii). Therefore, in this paper we concentrate 
on Ending the region of r and R in the simplified and tight situation. 

3.2 The result 

The procedure of our analysis is as follows: For a given r > 0 and a combination of 
signs of fermion masses, and M^r/R are calculated. We compute the neutrino mass 
matrix by the see-saw mechanism. Then, we compute neutrino mixings and neutrino 
mass ratios to see whether the result reproduce the desired ones. . 

We are fortunate that we found many regions of r which reproduce the desired mixings 
and masses. The result is summarized in Table 1 and 2. 
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In Table la, we show the result for mhrrir < 0. Solutions exist for various combinations 
of signs of quark and charged lepton as we see from the table. Each sign combination is 
expressed by the following abbreviated notation, Ei-Dj-UkX. The indices i, j and k runs 
1, 2, 3 and 4 which express sings of masses as (—, —), (—, +), (+, —) and (+, +) for a pair 
(e,/i) (denoted by E), (d, s) (denoted by D) and (m, c) (denoted by U), respectively, and 
X takes P and M that means the sign of mt is positive and negative, respectively. The 
region of r which reproduces the data in eqs.(23) is given in the column r. the column 
sin^ 2d^T- show the largest value of it, which is achieved by a specihc value of r shown in 
the column i?(10^^)(r). The value of R fixes the overall normalization of neutrino mass, 
especially niyr- We found solutions only for rrir < 0 in the range of r = 54 ~ 64 which 
coincide with our expectation given in Eq.(20). For all cases My£, has a non-hierarchical 
form, while MuR = RM takes a hierarchical form. The mixings and masses with specihc 
values of r and R for all cases are given in Table lb. 

In Table 2a, we show the result for mbfrir > 0 case. The notation to discriminate 
models which differ by sign combinations of fermion masses is the same as in Table la. 
We see that for all solutions nit < 0 and r ~ (500 ~ 960), except one case (E1-D2-U2P). 
These cases coincide what we expected from Eq.(21) and both Myj:, and M^r = RM take 
non(less)-hierarchical forms. 

The case E1-D2-U2P is realized when r ~ 35, which is also what we expected in 
Eq.(22), although the value of r is unknown. In this case, M^£, has a hierarchical form, 
while MyR = RM takes a non(less)-hierarchical form. In Table 2b, we show the mixings 
and masses for all cases. 

For completeness, the neutrino mass matrix and the neutrino mixing matrix for all 
cases are given in Appendix B, except for some specihc cases which we show below to see 
the details of the models. 

(1) nibnir < 0 case: 

(1-1) E1-D3-U2P with r = 58.2 
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The Dirac mass term and the Majorana mass term are given by 



/-0.0190 

-0.3530 

-0.1402 \ 

/ 0.0003 

0.0061 

0.0024 \ 

— 

-0.3530 

-2.2871 

-1.8954 , = R 

0.0061 

0.0445 

0.0326 


(-0.1402 

-1.8954 

2.6852 / 

(0.0024 

0.0326 

2.1704/ 


( 25 ) 

As we expected, has a non-hierarchical form the part relevant to the second and the 

third generations, while is hierarchical, because there is no cancellation of the 3-3 

element of The neutrino mass matrix and the neutrino mixing matrix is given by 

/1.165 20.46 8.402 \ / 0.964 0.067 -0.256 \ 

= -4 20.46 114.3 106.7 , = -0.242 0.613 -0.752 . (26) 

R 

V8.402 106.7 59.82/ V 0.107 0.787 0.608 

Then, we have sin^26*e^ = 0.018 and sin^ 26*^,- = 0.98. With R = 1.09 x 10^^, we found 

mi = -3.22 X IQ-^eV , m 2 = 2.25 x 10“W , m 3 = -1.83 x lO^^eV , (27) 

which are multiplied by 1.09 x 10^^/i? if the R dependence is kept. Thus 

Am ^2 = 5.08 X lO'^^eV^ , Am^g = 3.30 x lO'^eV^ . (28) 


As we have mentions before, the mixing angles are computed by assigning sin^e/i = 
{Uu)i 2 and sindf^T = {Ui,) 2 ^- This approximation is reasonable for the mixings and the 
hierarchical mass spectrum given in Eqs.(26) and (27). 

The sensitivity of mixing angles sm6fj_r and sinOefj. and also the ratio Ami 2 /Am |3 
with respect to the parameter r is shown in Figs, la, lb and Ic, respectively. The angle 
sin^^T- takes values larger than 0.7 for a wider range of r. The same holds for the ratio 
Ami2/Am23. The allowed region of r is essentially fixed by sin 6^611 and the region is 
between 58 to 62. This situation holds for all cases El-D3-UiP (i = 1 ~ 4). We will 
show the comparison with CHOOZ data, the disappearance test of z/g, at the end of this 
section. 

(1-2) E4-D2-U3P with r = 58.1 
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The Dirac mass term and the Majorana mass term are similar to the case (1-1), so 
that we give only the neutrino mass matrix and the neutrino mixing matrix 


/ 1.115 

19.83 

-9.592 \ 


/ 0.960 

0.0708 

0.270 \ 

19.83 

107.0 

-98.81 

,u,= 

-0.254 

0.620 

0.742 

V-9.592 

-98.81 

60.02 ) 


^-0.115 

-0.781 

0.614/ 


Then, we have sin^ 29e^ = 0.0200 and sin^ 26*^,- = 0.990. With R = 0.865 x 10^^, we found 
mi = -1.33 X IQ-^V , m 2 = -2.23 x lO'W , m 3 = 2.17 x lO-^eV , (30) 

which are multiplied by 0.865 x 10^^/i? if the R dependence is kept. Thus 

Am?2 = 4.97 X lO'^^eV^ , Am ^3 = 4.05 x lO'^eV^ . (31) 


As we can see from Figs.2a, 2b and 2c, the angle sin^^T- takes values larger than 0.7 
for a wider range of r. The same holds for the ratio Ami 2 /Am 23 . The allowed region of 
r is essentially fixed by sin^e/i- There are three allowed regions. This situation holds for 
all cases E2-D2-UiP, E2-D4-UiP, E4-D2-UiP, and E2-D4-UiP {i = 1 ~ 4). 

( 2 ) nibnir > 0 case: 

(2-1) E1-D3-U4M with r = 680 

The Dirac mass term and the Majorana mass term are given by 



/-0.2089 

-4.124 

-1.638\ 


/ 0.0003 

0.0061 

0.0024 \ 

— 

-4.124 

-29.95 

-22.15 

, — R 

0.0061 

0.0445 

0.0326 


^ -1.638 

-22.15 

-11.74 j 


U-0024 

0.0326 

-0.172/ 


As we expected, both M^j:) and M^ji have non-hierarchical forms for the part relevant to 
the second and the third generations. The neutrino mass matrix and the neutrino mixing 
matrix is given by 


/-141.3 

-2805 

-1111 \ 


/ 0.968 

0.0485 

-0.246 \ 

-2805 

-20159 

-15018 

,u,= 

-0.231 

0.553 

-0.800 

V -1111 

-15018 

-7336 j 


^ 0.0973 

0.831 

0.547 / 


( 33 ) 
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Then, we have sin^ O^fj, = 0.009 and sin^ 9^r = 0.92. With R = 1.22 x 10^^, we found 
mi = 2.21 X 10“W , m2 = 2.24 x IQ-^eV , m3 = -2.49 x lO'^eV . (34) 

Thus 

Amj^ = 5.00 X 10-®eV2 , Am^g = 16.15 x 10"W . (35) 


As we can see from Figs.3a, 3b and 3c, the angle sin6*^T- takes values larger than 0.7 
for a wider range of r. The same holds for the ratio Am^ 2 /^^ 23 - allowed region 
of r is essentially fixed by sin^e/i- There are three allowed regions similarly to the case 
(1-1). This situation holds for all cases El-D3-UiM (i = 1 ~ 4). 

(2-2) E3-D3-U3M with r = 790 

The Dirac mass term and the Majorana mass term are similar to the case (2-1) so 
that we only show the neutrino mass matrix and the neutrino mixing matrix 


214.7 

-3785 

-1502 \ 

/ 0.962 

0.0630 

-3785 

-28242 

-20425 

, = -0.257 

0.544 

-1502 

-20425 

-9864 j 

V 0.0946 

0.837 


As we can see from Eigs.4a, 4b and 4c, the angle sin^^T takes values larger than 0.7 
for a wider range of r. The same holds for the ratio Am^ 2 /^^ 23 - The allowed region of 
r is essentially fixed by sin 6^6^. There are two allowed regions. This situation holds for 
all cases El-Dl-UiM, E3-Dl-UiM, E3-D3-UiM (i = 1 ~ 4). 

(2-3) E1-D2-U2P with r = 36.5 

This is a very special case corresponding the case (b-2). The Dirac mass term and 
the Majorana mass term are given by 


MiyD 

MyR 


/-0.01399 0.21534 -0.10180 

0.21534 -3.1557 -1.1315 

V-0.10180 -1.1315 135.29 

/ 0.0003549 -0.0058998 0.002789 

R -0.0058998 0.094732 0.031000 

I 0.002789 0.031000 -0.17239 



(37) 
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As we expected, M^r, has a hierarchical form the part relevant to the second and the 
third generations, while is less hierarchical, because of the cancellation of the 33 


element in The neutrino mass matrix and the neutrino mixing matrix is given by 


/-0.559 

8.01 

28.0 \ 


( 0.999 

0.0374 

-0.00658 \ 

8.01 

-107 

-606 

, u,= 

-0.0290 

0.864 

-0.502 

1 28.0 

-606 

-800 J 


( -0.007 

0.502 

0.865 / 


Then, we have sin^ 6eu = 0.006 and sin^ Oar = 0.75. With R = 7.85 x 10^^, we found 


= -9.47 X lO'W , m 2 = 3.12 x lO'W , m 3 = -1.47 x lO'^eV . (39) 


Thus 


Am?2 = 9.71 X 10"®eV2 , Am^3 = 2.06 x lO^^eV^ . (40) 


As we can see from Figs.5a, 5b and 5c, sin is insensitive to r and the allowed region 
is determined by sin 6*^,- and Amf2/Am23. The allowed region of r is a tiny region, so 
that the solution is quite sensitive to the value of r. This case will be the most unlikely 
one. 

Finally, we comment on the CHOOZ data[2]. Since Aml 2 << Am23, consider 

the mixing angle for the CHOOZ disappearance test of may be dehned by sin6*ex = 
{Uu) 13 which is singer in our cases and the mass squared difference is AnR ~ = m|. 

Then, we compared these values of parameters with the CHOOZ two neutrino analysis. 
We found the followings: All cases except E1-D2-U2P predict sin26'er rsj ^0.2 0.3). For 

mbUir < 0 with m^rric < 0 cases and also mbnir > 0 cases, we predict Am? < 10“^eV^, so 
that they are safe. For mbm^ < 0 with m^mc > 0 cases, we predict Am^ ~ 2 x 10“^eV^ 
and their values are on the boundary of the excluded region. This can be remedied by 
changing r slightly from the value we took to achieve the largest value of sin 26 or by 
taking larger value of R to reduce the overall neutrino mass scale. For the E1-D2-U2P 
case, we predict sin26*eT ~ 1 x 10“*^, so that this satishes the bound. In summary, many 
of our cases predict in general rather large values of (77y)i3 ~ 0.25. However, all cases 
satisfy the CHOOZ bound. In several cases corresponding to mbmr < 0 with m^mc > 0, 
the z/g — oscillation is large enough to be observed in the near future experiments. 
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We explicitly showed some of our solutions. Other solutions show the similar matrices 
to one of the above cases. We listed the neutrino mass matrix rrii, and the neutrino mixing 
matrix Uiy in Appendix B for some of other cases. 

4 The summary 

We showed the mechanism to induce the non-hierarchical neutrino mass matrix by using 
the hierarchical forms of mass matrices of Md and Mi, although it is a quite non 
trivial problem. Our model contains only one adjustable parameter r which determines 
neutrino mixings and the ratios of neutrino masses. The other parameter R is used to 
determine the absolute magnitude of neutrino masses, so that it determines rriiyr in our 
hierarchical mass spectrum of neutrinos. 

Following our mechanism, we found many solutions, which are classihed into three 
cases. Our mechanism guarantees the non-hierarchical structure only for the part relevant 
to the second and the third generations and thus we inevitably predict the small mixing 
between z/g and z/^, while the large mixing between z/^ and z/^. We examined the r 
dependence of mixing angles and the ratios of neutrino masses and show the sensitivity 
of these quantities to r. Depending on the choice of relative signs of fermion masses, the 
pattern of sensitivity changed. We showed the solutions are not very sensitive to r except 
the case E1-D2-U2P case. 

There arises a question whether the existence of the solutions depends crucially on 
the values of quark masses and quark mixings, i.e., the value of tan/?. We analyzed by 
taking another set of values and found that the solutions given in this paper exist by 
the small change of the values of r and R. Thus the existence of our solutions does not 
depend on them. 

Finally, we comment on the CP violation. In the present analysis, we ignored the 
CP violation effects. When the CP violation tnrns on, many phases enter in onr model. 
One is from the Kobayashi-Maskawa phase and others are from the phases of qnark and 
charged lepton masses. These phases will relax the tight sitnation which we considered. 


17 



Thus, we expect that with the CP violation models cover the broader range of the mixing 
angles and the neutrino masses. This problem is under considerations. 

Also, there is no good reason to assume that Mi is diagonal, although it would be a 
hierarchical form. If we relax this assumption, the models would cover wider range of 
mixings and masses than what we obtained in this paper. 
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Appendix A: The proof of the incompatibility of the 
minimal model and quark and charged lepton mass 
spectrum 

For one 10 dimensional Higgs field n = 1 case, we have additional relation, 

Ml = CuDu + CdMd , (A.l) 


where 



(A.2) 


In this section, we hrst discuss that the model contains one 10 and one 126 repre¬ 
sentation can not reproduce the charged lepton masses and thus the model is rejected 
independently of neutrino masses and mixings. Then, we go to the model which contains 
n (> 2) 10 and one 126 models. We consider the ordinary see-saw mechanism to derive 
the left-handed neutrino mass matrix. In the case of type II see-saw model proposed by 
Mohapatra and Senjanovic[10], there appears the left-handed neutrino mass matrix. In 
our case, we do not consider this by taking vl is small. 

This model was considered by Mohapatra and his coworkers[4],[5],[8] as a model to give 
an unihed model to explain quark and lepton masses and mixings. Recently, Brahmachari 
and Mohapatra[8] showed that this model can not explain the desired neutrino mixing 
pattern, the small mixing between z/g and z/^, and the large mixing between z/^ and z/^-. 

Here, we show that this model is not able to explain the charged lepton masses, so that 
the model is rejected independently of the arguments on neutrino masses and mixings. 
Firstly, we see that the charged lepton mass matrix is written in a linear combination of 
quark mass matrices as in Eq.(A.l). 

Because \mr\ — \mi,\ ~ X^\mt\ with the the Cabibbo angle A = 0.2205, Cu ~ A^ and 
Q ~ 1. Knowing that Du is a diagonal matrix and Md is almost diagonal, and also 
the mass hierarchy of u-quark sector is much severer that that of d-quark, we observe 
that the contribution to the hrst and the second generation part of Mi from the u-quark 
part Du is negligible so that it is proportional to that of M^. Thus, the mass matrix 
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that predicts nie/nifj, ~ md/nis does not reproduce the observed hierarchical structure 
of d-quark and the charged lepton masses such as predicted by the Georgi-Jarskog mass 
relations rub = mr, rris = m^/3 and = 3me at GUT scale. 


Appendix B: The expressions of neutrino mass matrix 
and the neutrino mixing matrix 

Here, we show the explicit forms of the neutrino mass matrix mu and the neutrino 
mixing matrix Uy for various cases. From these matrices, the mixing angles and the ratios 
of neutrino masses are obtained. The result is summarized in Table lb and 2b. 


B.l mhUir < 0 case 


The case E2-D2-U1P with r = 56.1 predicts 


/-1.234 

18.60 

-8.698 \ 


/ 0.957 

0.0686 

0.283 \ 

18.60 

100.4 

-87.78 

,Uu = 

-0.266 

0.599 

0.755 

V-8.698 

-87.78 

48.17 j 


^-0.118 

-0.796 

0.592/ 


The case E2-D4-U1P with r = 62.9 predicts 


/-11.55 

21.09 

-10.92 \ 


/ 0.963 

0.0601 

0.263 \ 

21.09 

131.2 

-110.3 

,Uu = 

-0.248 

0.578 

0.778 

V-10.92 

-110.3 

52.03 j 


iy-0.105 

-0.814 

0.571 / 


(B.l) 


(B.2) 


The case E4-D4-U1P with r = 55.3 predicts 


/-6.952 

16.31 

-8.429 \ 


/ 0.962 

0.0594 

0.266 \ 

16.31 

97.70 

-82.65 

,Uu = 

-0.250 

0.583 

0.773 

iy -8.429 

-82.65 

40.67 ) 


^-0.109 

-0.810 

0.576 / 


B.2 The mbUir > 0 case 

(2-1) E1-D1-U2M with r = 830 
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/ 1529 

-3784 

-1677 \ 


/ 0.965 

0.0634 

-0.254 \ 

-3784 

-30058 

-22307 

, u,= 

-0.247 

0.545 

-0.801 

V-1677 

-22307 

-10215 j 


^ 0.0877 

0.836 

0.542 / 


The case E3-D1-U2M with r = 840 predicts 


2025 

-3876 

-1717 \ 


/ 0.964 

0.0654 

-0.258 \ 

-3876 

-30793 

-22848 

,u,= 

-0.252 

0.541 

-0.802 

-1717 

-22848 

-9985 ) 


^ 0.0872 

0.838 

0.538 / 


(B^5) 
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Table la: The ranges of parameters of r and R to predict the desired values of neutrino 
mixing angles and masses for mbTUr < 0. 


El-D3-UiP (i=l,2,3,4) 


e, /i 

d, 5 

u, c, t 

r 

R (10i3) (r) 

sin^ 26ij,t 


(+,-) 


58.9-60.9 

0.62-0.73 (60.9) 

■SB 



(-,+),+ 

58.2-61.7 

0.78-1.4 (58.2) 




(+,-),+ 

59.3-60.5 

0.61-0.78 (60.5) 




(+,+),+ 

58.5-61.4 

0.83-1.4 (58.5) 



E2-D2-UiP (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (10i3) (r) 

sin^ 26ij,t 

(-,+) 

(-,+) 


55.0, 55.8-56.1, 61.9-62.3, 63.3-63.9 

0.58-1.05 (56.1) 




(-,+),+ 

55.0-55.6, 56.4-56.6, 61.4-61.7, 62.6-63.8 

0.62-0.73 (61.4) 




(+,-),+ 

55.0-55.1, 55.8-56.1, 61.8-62.2, 63.2-63.9 

0.58-1.05 (56.1) 




(+,+),+ 

55.0-55.7, 56.5-56.7, 61.3-61.6, 62.5-63.7 

0.62-0.69 (61.3) 



E2-D4-UiP (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (10i3) (r) 

sin^ 26f tr 

(-,+) 

{+,+) 


55.0-55.1, 62.9-63.1, 63.8-64.4 

0.73-1.32 (62.9) 




(-,+),+ 

55.0-55.2, 55.8, 62.3-62.4, 63.0-64.1 

0.64-0.81 (62.3) 

mmm 



(+,-),+ 

55.1-55.2, 62.9-63.1, 63.8-64.4 

0.54-0.97 (55.2) 




(+,+),+ 

55.0-55.3, 55.8-55.9, 62.2-62.3, 63.0-64.1 

0.64-0.75 (62.2) 



E4-D2-UiP (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (10i3) (r) 

sin^ 26^r 

(+,+) 

(-,+) 


55.3-55.6, 56.8-57.7, 60.2-61.2, 62.6-63.2 

0.65-1.17 (57.7) 




(-,+),+ 

55.4-56.2, 57.2-57.6, 60.4-60.8, 61.9-63.1 

0.38-0.68 (60.8) 

WBBM 



(+,-),+ 

55.3-55.8, 57.0-58.1, 59.8-61.0, 62.4-63.2 

0.62-1.11 (58.1) 




(+,+),+ 

55.4-56.4, 57.4-57.8, 60.2-60.7, 61.8-63.1 

0.36-0.64 (60.7) 



E4-D4-UiP (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (10i3) (r) 

sin^ 26^t 

(+,+) 

(+,+) 


55.2-55.3, 62.7-62.9, 63.6-64.3 

0.54-0.98 (55.3) 

■EB 



(-,+),+ 

55.0-55.4, 55.9-56.0, 62.0-62.2, 62.8-64.0 

0.64-0.72 (62.0) 




(+,-),+ 

55.2-55.3, 62.6-62.9, 63.6-64.3 

0.67-1.21 (62.6) 




(+,+),+ 

55.0-55.4, 56.0-56.1, 62.0-62.1, 62.8-64.0 

0.53-0.58 (56.1) 
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Table lb: Predictions of neutrino mixing angles and masses for mbfrir < 0. All cases 
correspond to Table la. 


Model 
(r, i?/ 10 l 3 ) 

sin^ 20e^ 

sin^ 20 ^t 

mi (10 ®eV) 

m 2 (10 ®eV) 

m 3(10 ^eV) 


Am| 3 (eV^) 

E1-D3-U1P 
(60.9, 0.675) 

0.019 

0.89 

-0.0250 

1.89 

-4.08 

3.56 X 10“® 

1.66 X 10“® 

E1-D3-U2P 
(58.2, 1.09) 

0.018 

0.98 

-0.0322 

2.25 

-1.83 

5.08 X 10“® 

3.30 X 10-"^ 

E1-D3-U3P 
(60.5, 0.695) 

0.020 

0.89 

-0.00946 

1.97 

-3.91 

3.87 X 10“® 

1.52 X 10“® 

E1-D3-U4P 
(58.5, 1.115) 

0.019 

0.98 

-0.0231 

2.35 

-1.80 

5.51 X 10“® 

3.19 X 10“"^ 

E2-D2-U1P 
(56.1, 0.815) 

0.019 

0.98 

-0.304 

-2.26 

2.06 

5.03 X 10-® 

4.21 X 10-4 

E2-D2-U2P 
(61.4, 0.675) 

0.014 

0.89 

-0.436 

-1.93 

4.09 

3.54 X 10“® 

1.67 X 10“® 

E2-D2-U3P 
(56.1, 0.815) 

0.015 

0.98 

-0.289 

-2.25 

2.06 

4.98 X 10“® 

4.21 X 10-4 

E2-D2-U4P 
(61.3, 0.655) 

0.017 

0.90 

-0.411 

-1.87 

4.21 

3.33 X 10“® 

1.77 X 10-® 

E2-D4-U1P 
(62.9, 1.025) 

0.014 

0.96 

-1.29 

-2.58 

2.06 

4.98 X 10“® 

4.18 X 10-4 

E2-D4-U2P 
(62.3, 0.725) 

0.010 

0.86 

-1.77 

-2.63 

3.92 

3.80 X 10“® 

1.53 X 10-® 

E2-D4-U3P 
(55.2, 0.755) 

0.018 

0.96 

-1.32 

-2.59 

2.08 

4.97 X 10“® 

4.28 X 10-4 

E2-D4-U4P 
(62.2, 0.695) 

0.019 

0.89 

-1.80 

-2.60 

4.09 

3.53 X 10-® 

1.66 X 10-® 

E4-D2-U1P 
(57.7, 0.91) 

0.020 

0.99 

-0.0287 

-2.23 

2.02 

4.96 X 10“® 

4.04 X 10-4 

E4-D2-U2P 
(60.8, 0.53) 

0.004 

0.85 

-0.128 

-2.23 

5.02 

4.96 X 10“® 

2.51 X 10-® 

E4-D2-U3P 
(58.1, 0.865) 

0.020 

0.99 

-0.0133 

-2.23 

2.17 

4.97 X 10“® 

4.65 X 10-4 

E4-D2-U4P 
(60.7, 0.50) 

0.003 

0.85 

-0.109 

-2.23 

5.30 

4.97 X 10“® 

2.81 X 10-® 

E4-D4-U1P 
(55.3, 0.76) 

0.014 

0.96 

-1.09 

-2.49 

2.09 

5.01 X 10“® 

4.29 X 10-4 

E4-D4-U2P 
(62.0, 0 . 68 ) 

0.019 

0.89 

-1.49 

-2.37 

4.16 

3.39 X 10“® 

1.72 X 10-® 

E4-D4-U3P 
(62.6, 0.94) 

0.020 

0.97 

-1.08 

-2.49 

2.24 

5.01 X 10“® 

4.96 X 10-4 

E4-D4-U4P 
(56.1, 0.555) 

0.018 

0.89 

-1.45 

-2.34 

4.26 

3.36 X 10“® 

1.81 X 10-® 
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Table 2a: The ranges of parameters r and R to predict the desired neutrino mixings and 
masses for mbnir > 0. 


El-Dl-UiM (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (lOi®) (r) 

sin^ 26^r 




510-520, 550, 830-840, 900-950 

0.55-0.99 (550) 

0.92 




510-520, 540, 830-840, 900-950 

1.30-2.36 (830) 

0.92 




510-520, 550, 830-840, 900-950 

0.55-0.99 (550) 

0.92 



(+,+),- 

510-520, 540, 830-840, 900-950 

1.30-2.36 (830) 

0.92 


El-D3-UiM (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (lO^®) (r) 

sin^ 20ij,r 


(+,-) 


550-580, 640-740, 840-910 

0.84-1.51 (690) 





560-580, 640-740, 840-900 

0.87-1.57 (680) 





550-580, 640-740, 840-910 

0.84-1.51 (690) 




(+,+),- 

560-580, 640-740, 840-900 

0.87-1.57 (680) 



E3-Dl-UiM (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (lOi®) (r) 

sin^ 29fj,T 

(+,-) 



510, 540, 840-850, 910-960 

1.36-2.47 (840) 

0.91 




510, 540, 840-850, 910-950 

1.38-2.51 (840) 

0.92 




510, 540, 840-850, 910-960 

1.36-2.47 (840) 

0.91 



(+,+),- 

510, 540, 840-850, 910-950 

1.38-2.51 (840) 

0.92 


E3-D3-UiM (i=l,2,3,4) 


e, fi 

d, s 

u, c, t 

r 

R (lO^®) (r) 

sin^ 26f tr 

(+,-) 

(+,-) 


520-540, 570-580, 790-810, 880-940 

1.12-2.04 (790) 





530, 570-580, 790-810, 880-930 

0.63-1.13 (580) 





520-540, 570-580, 790-810, 880-940 

1.12-2.04 (790) 




(+,+),- 

530, 570-580, 790-810, 880-930 

0.63-1.13 (580) 



E1-D2-U2P 


e, fi 

d, s 

u, c, t 

r 

R (lOi®) (r) 

sin^ 26ij,t 


(-,+) 


36.2-36.5 

7.74-7.96 (36.5) 

0.75 
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Table 2b: Predictions of neutrino mixing angles and masses for mbfrir > 0. All cases 
correspond to Table 2a. 


Model 
(r, R/lQl®) 

sin^ 

sin^ 2e^T- 

mi (10 ®eV) 

m2 (10 ®eV) 

m3(10 ^eV) 

Amj2 (eV^) 

Am^g (eV^) 

El-Dl-UIM 
(550, 0.77) 

0.018 

0.92 

0.976 

2.44 

-2.63 

4.99 X 10-® 

6.90 X 10-4 

E1-D1-U2M 
(830, 1.83) 

0.016 

0.92 

0.941 

2.43 

-2.45 

5.01 X 10“® 

5.96 X 10-4 

E1-D1-U3M 
(550, 0.77) 

0.019 

0.92 

0.979 

2.44 

-2.64 

4.99 X 10-® 

6.90 X 10-4 

E1-D1-U4M 
(830, 1.83) 

0.016 

0.92 

0.942 

2.43 

-2.45 

5.01 X 10“® 

5.96 X 10-4 

E1-D3-U1M 
(690, 1.175) 

0.009 

0.91 

2.30 X 10-"‘ 

2.24 

-2.72 

5.00 X 10-® 

7.33 X 10-4 

E1-D3-U2M 
(680, 1.22) 

0.009 

0.92 

2.34 X 10-"‘ 

2.23 

-2.49 

4.99 X 10“® 

6.15 X 10-4 

E1-D3-U3M 
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Figure captions: 


Fig.l: The parameter r dependence of neutrino mixing angles and the ratio neutrino 
masses squared differences for E1-D3-U2P case {rribmr < 0). Fig.la, lb and Ic 
show the r dependence of sin^ 26 '^t-, sin^ 29efj, and Am^ 2 /^^ 23 - The allowed region 
of r is essentially determined by reproducing sin^ 29e^. 

Fig.2: The same as Fig.l for E4-D2-U3P case {mbfrir < 0). The allowed region of r is 
essentially determined by reproducing sin^ 29e^. 

Fig.3: The same as Fig.l for E1-D3-U4M case > 0). The allowed region of r is 

essentially determined by reproducing sin^ 29e^. 

Fig.4: The same as Fig.l for E3-D3-U3M case {mi,mr > 0). The allowed region of r is 
essentially determined by reproducing sin^ 29 

Fig.5: The same as Fig.l for E1-D2-U2P case {mbfrir > 0). The allowed region of r 
is essentially determined by reproducing sin^ 26^^.^ and Amf 2 /Am 23 . The allowed 
region is a tiny one, so that we need a hne tuning of r. 


27 






























